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ABSTRACT 

The timescales to replenish dust from the cool, dense winds of Asymptotic 
Giant Branch stars are believed to be greater than the timescales for dust de- 
struction. In high redshift galaxies, this problem is further compounded as the 
stars take longer than the age of the Universe to evolve into the dust production 
stages. To explain these discrepancies, dust formation in supernovae (SNe) is 
required to be an important process but until very recently dust in supernova 
remnants has only been detected in very small quantities. We present the first 
submillimeter observations of cold dust in Kepler's supernova remnant (SNR) 
using SCUBA. A two component dust temperature model is required to fit the 
Spectral Energy Distribution (SED) with T^arm ~ 102K and Tcoid ~ 17K. The 
total mass of dust implied for Kepler is ~ IMq - 1000 times greater than previ- 
ous estimates. Thus SNe, or their progenitors may be important dust formation 
sites. 

Subject headings: ISM: abundances — dust, supernova: individual — Kepler 

1. Introduction 

Recent blank field submillimeter (submm) surveys and observations of distant quasars 
have discovered a population of extremely dusty objects at high redshifts (Small, Ivison & 
Blain 1997; Bertoldi et al. 2003) implying that large quantities of dust are present in the 
Universe at z > 4 (Eales et al. 2003; Ivison et al. 2003, Issak et al 2002). Dust is generally 
believed to form in the cool atmospheres of giant stars and expelled via slow dense winds, yet 
at these redshifts, the time needed for such a process is greater than the age of the Universe. 
Morgan & Edmunds (2003) constructed models of the maximum dust mass produced in high 
redshift galaxies. Using a simple method to estimate the condensation efficiency of dust in 
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Stellar Winds (SW's) they found there is not enough time at redshifts > 4 for significant 
amounts of dust to be made in SW's unless star formation rates are very high and dust 
destruction rates are low. Dust production in SNe however provides an explanation for the 
dust seen in high redshift SCUBA sources as SNe evolve on a much shorter timescale. 

Theoretical models of nucleation of dust grains in supernova gas suggest condensation 

efficiencies of ~ 0.1 (Clayton, Deneault & Meyer 2001) and ~ 0.2 — 0.8 (Todini & Ferrara 
2001) could be reached in Type 11 SNe regardless of initial metallicity. However, the obser- 
vational evidence that SNe produce large quantities of dust is very weak. Infra-Red (IR) 
observations of young SNRs with IRAS and ISO found only ~ 10~^ and 1O~^M0 of dust 
in Tycho and Kepler (Douvion et al. 2001b) and 10^"^ in Cas A (Douvion, Lagage & 
Pantin 2001a; Arendt, Dwek & Moseley 1999). Observations of SN1987A have found a sim- 
ilar amount of dust (Dwek et al. 1992; Wooden 1997) although this could be much higher 
if the ejecta is clumpy (Lucy et al. 1991). One possible reason for the discrepancy between 
observations and theory would be if the dust in the SNRs is cold and virtually impossible to 
detect at IR wavelengths but observable in the submm waveband. Recent SCUBA observa- 
tions of Cas A (Dunne et al. 2003) have shown that a colder population of grains docs exist 
with an estimated 2 — 4Mq of dust, which implies a condensation efficiency of 0.5 — 0.8 for 
a 30Mq progenitor star. In this Letter, we present the first submm observations of another 
young Galactic SNR, Kepler's SNR. 

Kepler's SNR is the remnant from the early 1600's (Green 2001). Recent estimates place 
it at a distance of 4.8 - 6.4kpc (Rcynoso & Goss 1999) so we adopt 5kpc. The progenitor 
and type of Kepler's SN are still hotly debated. It was originally thought to be a Type la 
SN due to its light curve, and its distance from the Galactic plane suggests the progenitor 
was an old population II object. However, slow moving optical knots present in the region 
around the SNR show substantial overabundance of nitrogen, which can be most naturally 
explained by the the GNO cycle in massive stars (Bandiera 1987; Borkowski, Blondin & 
Sarazin 1992). The slow expansion velocity and high densities of the knots suggest they 
were not ejected in the SN explosion but rather in an earlier stellar wind, which again 
suggests a massive star progenitor (Blair, Long & Vancura 1991). X-Ray models of the 
remnant provide further evidence for a dense medium. These require ambient densities of a 
few H atoms per cubic centimeter (Kinugasa & Tsunemi 1991), which is 100 times greater 
than the ambient density of the interstellar medium (ISM) expected at 600pc above the 
Galactic plane (Whittet 2001). All these observations indicate there was dense circumstellar 
material (GSM) predating the SNR, which is most easily explained by the stellar winds from 
a massive star. For these reasons, Kepler's supernova has been classified as a Type lb SN. 
Bandiera (1987) has proposed that the progenitor of Kepler's SNR was a massive star which 
was ejected from the galactic plane. In this model, the star moves through its own stellar 
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wind, compressing the wind material into a bow shock in the direction of motion, and this 
produces the N-S asymmetry seen in the X-Ray and radio observations. 

2. Data Reduction 

We observed Kepler's SNR with SCUBA at the 15m James Clerk Maxwell Telescope 
(JCMT) on Mauna Kea in Hawaii. SCUBA simultaneously observes at 450 and 850/xm 
using bolometer arrays of 91 and 37 bolometers respectively. The field of view (fov) is ~ 2.3 
arcmin, but is slightly smaller at 450//m. The beam size (FWHM) at these wavelengths is 
~ 8 and 15 arcsecs respectively. As Kepler's SNR is actually larger than the fov of SCUBA, 
we observed the remnant with six 'jiggle-map' observations, each jiggle map being centered 
around the SNR with a chop throw of 180 arcsec. 

The data was reduced using the SURF software package (Sandcll ct al. 2001). We made 
corrections for atmospheric absorption using the opacities derived from skydip measurements 
with the JCMT and from polynomial fits to the sky opacity measured by the Caltech Sub- 
millimeter Observatory at 225GHz. Noisy bolometers were flagged and large spikes were 
removed. Any noise correlated across the array (i.e sky noise) was removed using the SURF 
program REMSKY. A map was made from the six individual datasets using one-arcsec pix- 
els. We created noise maps at each wavelength using the Monte-Carlo technique described 
in Eales et al. (2000). The S/N images at 850 and 450//m are shown in Figure 1 with 
the 850/um map clearly showing a ring like morphology similar to the X-Ray and radio ob- 
servations. The 450//m map shows emission in both the north-western and south-eastern 
quadrant. The noise at 450^m has a much stronger dependence on the atmospheric condi- 
tions than the noise at 850;um. The conditions were much worse when we observed the other 
two quadrants, so the lack of obvious emission in those quadrants and also in the centre of 
the remnant may partly be the result of the increased noise in these regions. 

The image was calibrated using primary calibrators which were observed with similar 
chop throw to the observations of the SNR. The flux was measured in eight apertures, which 
were placed on the maps to include all the obvious emission from the SNR (using the S/N 
maps as a guide). The flux calibration was carried out using apertures of the same size on 
the calibrator map. The errors in the fluxes were estimated as described in Dunne et al. 
(2000) with calibration errors of 15% and 20% respectively. The fluxes measured over the 
whole remnant at 850 and 450/im are 1.0 ± 0.2Jy and 3.0 ± 0.7Jy respectively. Thus the 
ratio of 450 to 850yum flux for the SNR as a whole is 3. 
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3. Estimating the Dust Mass 

Supernova remnants are strong sources of synchrotron radiation which can contaminate 
the longer submm wavelengths. To obtain the submm flux which is emitted by dust grains 
only, we need to correct for this contamination. Using the average of the IR fluxes from the 
literature (Arendt 1989; Saken, Fesen & ShuU 1992; Braun 1997) and the submm fluxes from 
this work, we can produce Kepler's IR- Radio SED (Figure 2). This clearly shows that the 
emission from dust is in excess of the synchrotron emission. We can estimate the contribution 
of the synchrotron to the submm using the flve radio measurements which arc well fitted by 
a power law slope S,^ cx z/" where a = —0.71. Hence the synchrotron emission is responsible 
for 30% of the submm flux at 850//m and 6% at 450//m. 

Comparison of the 450 and 850//m maps with a radio map also clearly indicates that 
there is submillimetre emission in excess of the expected synchrotron emission. Figure 1(c) 
shows the 5GHz VLA image of the SNR (Delaney, private communication). The map shows 
the well known N-S asymmetry associated with the interaction with a dense CSM. The 
850/im map also shows a N-S asymmetry but is much less pronounced than the radio. The 
450//m map shows an even lesser N-S asymmetry. A point-by-point comparison of the maps 
also indicates excess submm emission. The strong emission visible in the south west on the 
submm maps is at a position where there is hardly any emission on the radio map. 

To correct the submm images for the contribution of synchrotron radiation we smoothed 
the VLA radio image to the resolution of the submm images and scaled it to the submm 
wavelengths using the power law spectrum shown in Figure 2. We then subtracted it from 
the submm image. The synchrotron-subtracted 850//m image is shown in Figure 1(d). The 
brighter northern limb in the 850//m unsubtracted image has been removed and the 450 
and 850/im images now follow a similar morphology with two bright emission regions in the 
north and south of the remnant. The peaks on the 450 and 850/xm maps do not exactly 
coincide with each other but this is not unexpected given the low S/N at 450yum and the 
fact that we arc extrapolating a long way in frequency from 5GHz to the submm waveband. 
The remnant is known to have variations in the spectral index of order ±15% (Delaney et 
al. 2002). This may produce errors in the synchrotron-corrected maps although the general 
agreement between the 450 and 850//m images suggests this is not a problem. Note also that 
the fact that the ratio of 450 to 850/xm is 3 makes it practically impossible to explain the 
excess submm emission by some synchrotron process. 

The dust mass can be measured directly from the flux at submm wavelengths using 
Mfi = ^ B^^uT) (Hildebrand, 1983) where S„ is the flux density measured at frequency u, 
D is the distance and Ki, is the dust mass absorption coefficient. B is the Planck function 
and T is the dust temperature. We tried fitting a single temperature greybody to the SED 
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{S^ oc B^ul^) but we could not obtain an acceptable fit. So we fitted a two-temperature 
grcybody allowing P and the warm and cold temperatures to vary. The best fit SED gave 
P = 1.2, = 102K and Tc = 17K. The possible reasons why the dust could be so cold are 
discussed in Dunne et al (2003) and will be considered in future work (Dunne et al, in prep). 
We cannot rule out the presence of dust existing at intermediate temperatures between these 
two extremes until further observations around the peak of the SED (200yLtm) are made. We 
used a bootstrap technique to derive errors on these values. We created 3000 sets of artificial 
fluxes from the original data and then fit our two temperature model to each set. We derived 
our errors from the distribution of T^, Tc and /? produced by these fits. 

The largest uncertainty in the dust mass comes from the uncertainty in k. We have 
followed Dunne et al. (2003) in trying three different values of k from the literature: (1) 
1^850 — 0.6 — l.lm^kg'^, the range observed in laboratory studies of clumpy aggregates. (2) 
^^850 = 0.16 — 0.8m^%~^ the range observed in circumstellar environments/reflection nebulae 
and (3) — 0.04 — 0.15m^%~^, the range observed for the diffuse ISM where dust is likely 
to have encountered extensive processing. In each case we used the average values of the 
ranges to estimate the dust mass which are given in Table 1. In Cas A, the higher k values 
were required to give a reasonable dust mass. Using the highest k values for Kepler gives 
a lower range of 0.2 — O.3M0 whereas using the lowest k values gives 1 — 3Mq. Whichever 
value of K we use we obtain at least three orders of magnitude more dust than obtained from 
mid-IR measurements. This is an order of magnitude lower than the dust observed in Cas 
A but the far greater brightness of Cas A than Kepler at all wavelengths imphes a (as yet 
not understood) difference in the dynamics of the explosions. 



4. Discussion 

Kepler's SNR still appears to be in an early evolutionary phase in which the swept 
up mass is not much more than the total ejected mass (Dickel 1987; Hughes 1999). If the 
interstellar density at the location of Kepler's SNR is uh ~ 0.001 — O.lcm"^ (Bandiera 1987; 
Borkowski et al. 1992) then given the size of the remnant and assuming a normal gas to 
dust ratio of 160, the maximum dust mass swept up by the SNR is ~ 10~^Mq. 

Evidence for an even denser CSM comes from slow moving optical knots observed out- 
side of the ejecta. These are thought to originate from accelerated clumps of circumstellar 
material and have densities Uh ~ lOOOcm"^, which imply a pre-shock density of lOOcm"^ 
(Blair et al. 1991). If this density is representative of the CSM surrounding the entire rem- 
nant, then the blast wave could have swept up IOOM0 of gas and a dust mass of ~ IMq. 
This is however, at odds with the mass observed in the X-Ray which suggests that there is 
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only ~ 2 — 5Mq of gas within the SNR (Kinugasa & Tsunemi 1991; Borkowski et al. 1992). 

Some additional arguments that the dust has not been swept up come from considering 
the morphologies of the different images. Some of the dust is found between the forward 
and reverse shocks of the X-Ray gas, where the gas is thought to be dominated by the 
ejecta rather than the swept up material (Hughes 1999). A large fraction of the dust is 

also found in the south of the SNR, where the radio image shows no evidence for a strong 
interaction between the SNR and a dense CSM. However, some of the submm emission is 
found outside the X-Ray and Radio remnants. This may be due to some of the dust being 
formed in the progenitor's stellar wind or due to the velocity of the ejecta moving faster in 
these regions. The lack of a detailed position for the two shock fronts does not allow for a 
complete comparison though further work will try to address this problem (Dunne et al, in 
prep). 

Could the dust be from a massive star, formed in its pre-supernova RSG and WR phases? 
The maximum observed dust mass lost in WR winds is 6 x 10~'^MQyr~^ (Marchenko et al. 
2002). The violent mass loss from the WR will blow the material outwards into a shell with 
outer radius depending on the velocity and duration of the wind. If the wind velocity is 
~ 1000A;ms~^ and if the WR wind was responsible for the dust seen at the radius of the 
submm emission, the WR phase of the precursor must have only lasted around 10^ years. 
This would produce only O.OIM0 of dust. The interaction between the previous RSG and 
WR winds may have created more dust but the maximum yield of heavy elements ejected 
from the most massive stars during the RSG phase is < 2Mq. The condensation efficiency 
in these atmospheres was estimated to be < 10% (Morgan & Edmunds 2003) and hence it is 
difficult to produce all of the observed dust in this way. Note, if the condensation efficiency 
is indeed higher, this could explain the lower end of the cold dust mass in Kepler and hence 
SNe progenitors would be responsible for polluting the high redshift galaxies. 

The mass of dust we observe is lower but similar to the mass of dust found in Gas A. If 
Kepler's progenitor was a SOM© star (Bandiera 1987), the condensation efficiency could have 
been ~ 0.1 — 0.6. If we assume that the dust in Kepler's SNR is representative of the dust 
formed in every 30 Mq progenitor then this is roughly comparable with the dust mass/SNe 
given in Todini & Ferrara (2001). Using their models with the condensation efficiency ob- 
tained for Kepler's SNR from this work, we estimate that SNe inject ~ 5 — 15(x lO^'^)M0|/r^^ 
into the ISM. If the dust injection rate for giant stars is ~ 5(xlO^"^)M0yr~^ (Whittet 2001; 
Morgan & Edmunds 2003) then SNe are injecting around 1-3 times more dust than stellar 
winds. Adding supernova or massive stars as a source of dust to the Galactic interstellar 
medium explains the high redshift submm observations and could well resolve the current 
discreprancy (Jones et al. 1994) between dust production and destruction rates. 
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Table 1: The parameters from the SED fitting using minimum fit to the points. The dust 
masses are the median of the distribution from the bootstrap technique (same as the best fit 
values) with 68% confidence interval quoted as the error. The values of k are (as in Dunne 
et al. 2003) Ki(450) = 1.5m'^kg~\ ki(850) = 0.76m'^kg-^ (laboratory studies); K2(450) = 
O.SSm^kg"^, fi;2(850) = O.Srri^kg^^ (observations of evolved stars); ft;3(450) = 0.26m^kg^^, 
Ks{850) = O.OTm^kg-^ (observations of the diffuse ISM). 
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Fig. 1. — Signal/Noise maps of the Kepler SNR. (a) 850yum smoothed to 14" with contours 
3a", 4(j, 5a"... etc. (b) 450/um image smoothed to the 850/um resolution with contours 1.5a", 
2.5a-, 3.5a-... (c) 6cm radio image of Kepler (Delaney et al. 2002) smoothed to the resolution 
of the 850//m beam, (d) The synchrotron subtracted 850/xm image with levels showing 3a-, 
4a-, 5a-... 
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Fig. 2.— Kepler's IR - Radio SED for the IRAS 12, 25, 60 and 100/xm fluxes, the 850 and 
450/xm measurements from this work and the radio (Matsui et al. 1984; Delaney et al. 2002). 
The solid hne represents a two temperature fit to the data points with P — 1.2, Ty^ — 102K, 
Tc = 17 K. The straight line through the radio measurements is a power law with a spectral 
index -0.71. The lOOyum flux was chosen here to be 5Jy rather than the lower value of ~3Jy 
observed by Arendt (1989) and Braun (1987). This was a conservative choice as using the 
lower value would produced a more pronounced second peak in the SED. To obtain the best 
fit greybody through these points would then require a cold temperature of ~ lli^ and 
hence more than double the derived dust mass. 



